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Abstract 

In this paper we study the Leray weak solutions of the incom- 
pressible Navier Stokes equation in an exterior domain. We describe, 
in particular, an hyperbolic version of the so called artificial compress- 
ibility method investigated by J. L. Lions and Temam. The convergence 
of these type of approximation show in general a lack of strong con- 
vergence due to the presence of acoustic waves. In this paper we face 
this difficulty by taking care of the dispersive nature of these waves 
by means of the Strichartz estimates or waves equations satisfied by 
the pressure. We actually decompose the pressure in different acoustic 
components, each one of them satisfies a specific initial boundary value 
problem. The strong convergence analysis of the velocity field will be 
achieved by using the associated Leray-Hodge decomposition. 
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1 Introduction 



This paper is concerned with the description of a type of approximation 
method for the weak solutions of the 3-D Navier Stokes equations over 
an exterior domain Q. We say that f2 is an exterior domain if it is the 
complement in M 3 of a compact set (usually called compact obstacle). In 
order to solve the Navier Stokes equations in an exterior domain we need to 
find the velocity and the pressure fields which together solve the equations 
and which moreover assume given boundary data on the obstacle. This 
latter requirement contains usually the most relevant difficulty in this type 



1 



of analysis. 

The most simple examples of such kind of flow are fluids filling up space 
and flowing past spheres, plates and cylinders as the flow of a river around 
the stones lying on the riverbed. The interest in studying the Navier Stokes 
equations in such kind of domains arises from many phenomena in physics 
and applications in engineering models (see [48J and references therein). 
More complicated examples are the motion of bubbles in a liquid ( we may 
think at the bubbles in the ocean) and the sedimentation of particles. In 
these cases it is important to determine the forces that the fluid exerts on 
the structures. Similar phenomena comes also from climate modeling as for 
the rain drops falling within clouds in the high atmosphere and for these 
reasons even the flow around the sphere, as the simples example of a falling 
particle, is still a subject of great interest [37]. Other examples can be 
included in the same framework from engineering models like the design of 
project aircrafts wings with a high speed airfoil. Interesting information 
concerning these questions can also be obtained by studying insect flight, 
|47| . Indeed, if we consider the air as a fluid the insect can be seen as a 
moving obstacle. These kind of measurements are considered important and 
have many applications in construction micro-air vehicles for reconnaissance 
mission and space missions as the flight in the atmosphere of Mars, A 
final example comes from the area of medicine in the study of hemodynamics 
where it could be of great importance to understand and model the blood 
flow around an embolus or a prothesis. 

The mathematical model of an incompressible fluid in an exterior domain is 
given by the following set of equations 



dtu + div(-u (g> u) — fiAu = Vp + f x £ ft, t > 

divu = x £ $7, t > 0, 

u(x, 0) = uq(x) x E fi, 
u\dn = 0, liniM^oo u = 0, 



(1.1) 



where Q, is an exterior domain of M 3 , u G R 3 denotes the velocity vector 
field , p £ M. the pressure of the fluid , / € K 3 is a given external force, 
\i is the kinematic viscosity. In the case of the whole space M. d or of a 
bounded domain there exists, in the mathematical literature, several results 
concerning the existence and regularity of Leray weak solutions to the Navier 
Stokes equations, for example we can refer to books of P.L. Lions [32] and 
Temam [46] , The exterior problem for the Navier Stokes equation consists 
of finding in the region exterior to a closed bounded surface, velocity and 
pressure functions which together solve the equations and are such that the 
velocity assumes given values on the surface and tend to a prescribed limit 
at infinity. For what concerns the existence theory of weak solutions in the 
case of an exterior domain a first result can be found in the paper of Leray 
[30j . We recall here what we mean for a Leray weak solution of the system 
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Definition 1.1. We say that u is a Leray weak solution of the Navier Stokes 
equation if it satisfies (jl.lj) in the sense of distributions, namely 

J J ^Vu • Vv — UiUjdnfj — u ■ dxdt 

(/, <p) H -i xH idxdt + / u ■ ipdx, 
o Jn 

for all ip £ C$°(n x [0,T]), div^ = and 

divu = in V(Q x [0,T]) 
and the following energy inequality holds 



u(x, t)\ 2 dx + I I | Vu(x, t)\ 2 dxds 
q Jo Jn 



<- / |uo| 2 (ix+ J (f,u) H -i xH ids, for all t > 0. 
2 Jo 

Similarly to the problems on the whole space or on bounded domains, 
also here, there is a large amount of literature concerning the regularity 
of solutions. For instance, Finn in [T3], [T3], [E] studied the exterior 
stationary problem within the class of solutions, called by him as "physically 
reasonable", which tend to a limit at infinity like |x| - 2 -? ? for some ij > 0. 
For small data he proved both existence and uniqueness in this class. In 
the case of nonstationary flow the question of existence and uniqueness was 
later addressed by Hey wood. In [23] he recovered stationary solutions as 
the limit of nonstationary solutions. Then, he also studied the stability 
and regularity of nonstationary solutions in |25j . [24] . [26] . [27] . Further 
regularity properties of solutions in exterior domain can be found in |18j . 
|19j . [33], [36], [29]. In [20] the theory of the Navier Stokes equations was 
developed considering a moving or a rotating obstacle. 

Motivated by the previous examples and applications there have been 
also considerable efforts to develop numerical approximation methods. One 
of the major difficulty regards the development of numerical schemes in- 
cluding in an efficient way the incompressibility constraints. Chorin [5], [6], 
Temam [44j . [45] and Oskolkov [38], in the case of a bounded domain, to 
overcome the computational difficulties connected with the incompressibility 
constraints introduced what they named "artificial compressibility approx- 
imation". They considered a family of perturbed systems, depending on 
a positive parameter e, which approximate in the limit the Navier Stokes 
equation and which contain a sort of "linearized" compressibility condition, 
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namely the following system 



< 



d t u £ + Vp £ = /iAu' 
edtp £ + divu e = 0, 




(1.2) 



where tt G R 3 , (x,t) G Q x [0,T], u e = n e (x,t) G M 3 and p £ = p £ (x,t) G R, 
f = f(M)e« 3 



The papers of Temam [33], [35] and his book [36] discuss the convergence 
of these approximations on bounded domains by exploiting the classical 
Sobolev compactness embedding and they recover compactness in time by 
the well known J.L. Lions [31] method of fractional derivatives. The same 
system was used in [TU] in the case of the Navier Stokes equations in the 
whole space R 3 and modified in a suitable way in [8] for the Navier Stokes 
Fourier system in R 3 . In those papers the authors have carefully to estimate 
the acoustic waves for the pressure p £ which is the cause of the lost of the 
strong convergence. In order to overcome these difficulties they exploit the 
dispersive properties of these waves, by using L p — L q estimates of Strichartz 
type [2T] , [28] . [4*3] . Here, in order to approximate the system (jl.ip we in- 
troduce the system (11.21) with appropriate boundary data, which makes the 
problem considerably more difficult as we will see later on. We can see from 
the second equation of (|1.2p that as e goes to 0, the acoustic pressure waves 
propagate with high speed of order 1/e in the space domain. Because of the 
fast propagating of the acoustics one expect the velocity u £ to converges only 
weakly to the incompressible solution of the Navier Stokes equation. Here 
we overcome this trouble by using the dispersion of these waves at infinity 
obtaining the strong convergence of u £ . In this paper we need dispersive es- 
timates of Strichartz type on exterior domains. These estimates have been 
recovered by Smith and Sogge, [31], [32] in the case of odd space dimension 
and by Burq [3] and Metcalfe [33] for even space dimension (see Section [2T2]) . 
The connection with the dispersive analysis of the acoustic wave equation 
has also been considered to study the incompressible limit problem. Similar 
phenomena appear also in the modeling the Debye screening effect for semi- 
conductor devices, [11] . It is worth to mention here that this type of singular 
limits from hyperbolic to parabolic systems is not covered and doesn't fits 
in the general framework of diffusive limits analyzed in [9] . 

In order to understand the additional difficulty with respect to the whole 
domain case it is important to remark that the presence of a boundary re- 
quires consistent boundary conditions for the approximating problems. This 
problem is not present for the approximating velocity u £ since the limit ve- 
locity satisfies Dirichlet boundary condition, hence we can impose the same 
condition on u £ . Unfortunately for the pressure we have a different situa- 
tion, indeed for the limit problem (incompressible Navier Stokes) the natural 
physical boundary condition is of Neumann type, while the wave equation 
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structure for p £ and the Strichartz estimates require Dirichlet type condi- 
tions. This problem is well known also in the numerical literature (see for 
instance Chorin [I], Gresho and Sani [22j . Sani, Shen and al. |39j and the 
references therein). Motivated by the formal analysis of the previous men- 
tioned papers we introduce here appropriate Dirichlet boundary conditions 
for the pressure which are expected to be consistent with the Neumann 
boundary condition of the limit problem (see Section 3 for further details). 
However, to make our analysis rigorous it will be necessary to decompose 
the approximating pressure in different acoustic components, each one of 
them with its own appropriate boundary conditions. 

This paper is organized as follows. In Section 2 we recall the mathe- 
matical tools and basic definitions that we need through the paper and we 
describe the Strichartz estimate we are going to use. In Section 3 we intro- 
duce the approximating system and we state our main result. The Section 
4 is devoted to the a priori estimates that are derived from standard energy 
type estimates. In Section 5 we recover further estimates by exploiting the 
wave equation structure for the pressure. In Section 6 we show the strong 
convergence of the approximating sequences. Finally, in Section 7 we prove 
our main result. 

2 Notations and Preliminaries 

For convenience of the reader we establish some notations and recall some 
basic theorems that will be useful in the sequel. 

From now on $7 denotes an exterior domain to a compact obstacle in M. d . 
Precisely, Q is the complement in R d to a compact, strictly convex, smooth 
set contained in {|x| < R}. Moreover, 0, is assumed to be non trapping in 
the sense that there is a number Lr such that no geodesic of length Lr is 
completely contained in {|x| < R} n fi. 

2.1 Definition of spaces and Operators 

We will denote by V{VL x M + ) the space of test function C^°(0 x R+), 
by P'(Q x M_|_) the space of Schwartz distributions and (•, •) the duality 
bracket between V and V. Moreover W k ' p (n) = (I - A)~£ L P ({1) and 
H k (Q) = W k,2 (Q) denote the nonhomogeneous Sobolev spaces for any 
1 < p < co and k E M. The notations L P L% and L P Wx ,q will abbre- 
viate respectively the spaces L p ([0, T];L q (Q.)) and L p ([0, T]; W k > q (Sl)). In 
the definition of homogenous Sobolev spaces we have to be more precise. 
The homogeneous Sobolev norm H^iW^ 1 ) on the whole space W 1 is given by 
ll/llifT (R d ) = ll("V^^) 7 /ll-L 2 (R d )- On the exterior domain Q we have to define 
the space £[^(£1), 7 £ R which is the homogeneous Sobolev space associated 
to the square root of the Laplace operator with Dirichlet boundary condition 
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on O, \J— Ad- To be more precise fix (3 G C^°(M. d ) a smooth cutoff function 
such that (3(x) = 1 for \x\ < R and let be Q a compact manifold with the 
boundary containing Br = Q n {\x\ < R}, we are able to define 

\\f\\Hi(n ) = \\m HU n) + \\( 1 -mH^y 
Notice that for functions with support in {|x| < R} we have 

H/llffJCn) = 11/11^2,(0)- 

Functions / G H^tl) satisfy the Dirichlet conditions f\g^ = and when 
7 > 2 we must require the compatibility condition 

Aj f\an = for 2 J < 7- 

With the Dirichlet condition fixed we may define the spaces if 7 ($7) in terms 
of eigenfunctions of A. Since Q is compact we have {vj} C H^O,) n C°°(fi) 
an orthonormal basis of L 2 (Q) with Avj = —XjVj, Xj > 0, Xj f oo. For 
7 > we define 

= {/ G L 2 (f2) | £ |/(j)| 2 Aj < +oo|, 
L i>o J 

where /(j) = (f,Vj). The f/^, norm is given by 

ll/ll| 7( n)=El/(i)l 2A J- 
i>o 

For 7 < 0, we define H^Q) in term of duality. Moreover we mention that 
for r < s, 

Wf\\H r D {n) - C 'II/II%(Q)- 
For further details see [3] and [M] . 



Let us define in the case of an exterior domain the Leray's projector P 
on the space of divergence - free vector fields and Q on the space of gradients 
vector fields. Any vector field F on Q with F\qq = can be decomposed 
in the form 

F = K + VA, 

where div K = and K ■ n = on dQ. (n is the exterior normal vector to 
In fact we can recover A as a solution of the following system 

'AA = divF 

— = in dQ, / A = 0, 
on J n 
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that is well known has a unique solution (see [7]). It turns out that P and 
Q assume the form 

Q = VA^div P = I-Q, (2.1) 

where A^ 1 is the inverse of the Laplace operator with Neumann boundary 
condition. It can be proved that P and Q are projections. Let us remark 
that contrary on what happens in the case of the whole domain, here P 
and Q don't commute with translations and so with derivatives, but they 
are still bounded in LP and so in W k ' p (1 < p < oo) space (see |12] and 
references therein). 

2.2 Strichartz estimates in exterior domain 

Let us consider the following wave equation defined in the space [0, T] x R d 

(<9 4 2 - A) w(t, x) =F(t,x) 
w(0,-) = f, d t w(0,-)=g, 

for some data f,g,F and time < T < oo. As is well known the wave 
equation belongs to the so called dispersive equations. In 1977 Strichartz 
|43| realized that combining the dispersive properties of the wave equation 
with the restriction theorem of the Fourier transform on manifolds he could 
set up the following estimate 

IML*!,* + W d t W \\ L 4 W -^ £ ll/ll + \\9\\h-1/2 + 11-^11^/3^4/3 • (2-2) 

Later on this estimate was generalized to the following one (see |21j . |28| ) 

\W\l<lz + WdtwW^-^r < \\f\\ m + \\g\\ m -i + \\n L f L r>, (2-3) 
where (q, r), (q, f) have to be wave admissible pairs, namely they satisfy 
2 , N /l 1\ 2 , . (\ 1\ 

1 d d Id 
- + - = 7^-7 = - + -- 2. 
q r 2 q' r' 

The estimate (|2.3|) still go under the name of Strichartz estimate. A further 
generalization of the estimate (12.30 is given when we consider the wave 
equation on an exterior domain with Dirichlet boundary conditions, namely 
w is a solution of the following system 

' (pi - A) w(t, x) = F(t, x) (t,x)£R + xQ 
i w(0,-) = f(x)eHl 
* d t w(0,x)=g(x)£Hl~ 1 
w(t, x) = 0, x £ dQ, 
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for some data f,g,F and time < T < oo. Then, w satifies the following 
Strichartz estimate, 

\\M\lili + \\ 9 M\ L?w -i,r < WfWfn + bllM- 1 + ll F IL«' L f" ( 2 - 5 ) 

provided that (q,r), (g, f) are wave admissible pairs in the sense of (I2.4p . 
As we can observe the estimate (|2.5p has the same structure as the one (|2.3|) 
on the whole space M d , but in order to prove it is necessary to establish new 
decay estimates. These estimates are different if the space dimension is odd 
or even because of the lack of strong Huygen's principle in the latter case. 
The local Strichartz estimate for the homogenous case was proved by Smith 
and Sogge in |41j . Then in [i2j they established for the nonhomogenous 
wave equation the global estimate in space and time for odd space dimen- 
sion. The even space dimension estimate was obtained independently by 
Metcalfe [M] and Burq [3]. 

Of course the estimate (|2.5|) includes the case of the original estimate ob- 
tained by Strichartz in 1977 [43] . namely 

IMIl*L< + \\ d M\ L 4 W -^ < H/IUl/2 + HslL-l/2 + ||F|| 4/3 4/3, (2.6) 
t x yy x D D t x 

From the estimate (|2.6p can be deduced the following one that will be usefull 
for us in the sequel, 

WMklLi + ll^ll^-M < WfW^, + ||5||^-i/ 2 + \\F\\l\lI- (2-7) 

Later on we shall also use (|2.5|) in the case of d = 3, (<f , f') = (1, 3/2), then 
7 = 1/2 and (q,r) = (4,4), namely the following estimate 

\\M\lI x + \\dtw\\ L , w -iA < ||/||^ /2 + \\g\\ 6 -i/> + \\n L i L l/*- (2-8) 
2.3 Preliminary theorems and lemma 

Finally we mention here the following technical lemma that will allow us to 
get LP estimate by means of W~ k)P norms. 

Lemma 2.1. Let us consider a smoothing kernel j £ Cq°(£1), such that 
j > 0, J^ d jdx = 1, and define the Friedrichs mollifiers as 

j a (x) = a~ d j Q . 
Then for any f £ H 1 ^), one has 

\\f - f * 3 a \\ L p(n) < CX^IIV/IU^Q), (2-9) 

where 

£>G[2,oo) ifd = 2, p£ [2,6] if d = 3 and a = d 
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Moreover the following Young type inequality holds 



\\f*ja\\ L nn)<Ca s d U *)||/|| w -.„ (n) , (2.10) 
for any p, q € [1, oo], q < p, s > 0, a G (0, 1). 

Finally, we need to recall the following compactness tool (see |40j). 



Theorem 2.2. Let be T C L p ([0, T]; B), 1 < p < oo, B a Banach space. T 
is relatively compact in L p ([0,T]; B) for 1 < p < oo, or in C([0, T];B) for 
p = oo if and only if 

(i) <| J f(t)dt, f 6 B j is relatively compact in B, < t\ < t2 < T, 

(ii) lim H/ (t + h) - f(t)\\ L pr [0T _ h] . B) = uniformly for any / Gf. 

h— >0 

3 Approximating system and main result 

As we explained in the Introduction in order to approximate the system 
we introduce the following system 

d t u £ + Vp £ = fiAu £ - (u £ -V)u £ - i(div u £ )u £ + f £ , 
edtp 6 + div u £ = 0, 

where (x,t) £ fix [0,T], u £ = u £ (x,t) £ R 3 and p £ = p £ (x,t) £ R, f £ = 
f £ (x,t) £ M 3 . As we can notice the constraint "divu = 0" of the system 
(jl.ip has been replaced by the evolution equation 

dtp £ = — div u £ , 
e 

which can be seen as the linearization around a constant state of the con- 
tinuity equation in the case of a compressible fluid. Concerning the first 
equation of the system (|3.1|) we can observe that compared to the equation 
of the balance of momentum it has the extra term — 1 /2(div u £ )u £ which has 
been added as a correction to avoid the paradox of increasing the kinetic 
energy along the motion. 

Since it will not affect our approximation process, for semplicity, from 
now on, we will take \x = 1 and f £ = 0. 

Furthermore we assign to the system (|3.ip the following two initial con- 
ditions 

u £ {x, 0) = u £ (x), p £ {x, 0) = p £ {x). 

It is worth to mention here that the Navier Stokes equations require only 
one initial condition on the velocity u. Hence our approximation will be 
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consistent if the initial datum on the pressure will be eliminated by an 
"initial layer" phenomenon which will be a consequence of the dispersive 
nature of the acoustic pressure waves. Since in the limit we have to deal 
with Leray solutions it is reasonable to require the finite energy constraint 
to be satisfied by the approximating sequences (u e ,p £ ). So we can deduce a 
natural behaviour to be imposed on the initial data (uq,Pq), namely 

u £ = u £ (-,0) — ► u = u(; 0) strongly in L 2 (Q) 
VePo = V£P e (-,0) — ► strongly in L 2 (Vt). 

Let us remark that the convergence of y/spQ to is necessary to avoid the 
presence of concentrations of energy. 

Since we are in an exterior domain the system (|3.ip needs to be supple- 
mented with boundary data. Taking into account the Navier Stokes equa- 
tions p.ip the natural choice is to assign homogenous Dirichlet boundary 
condition to the velocity vector field u £ , namely 

u £ \ dn = 0. (BC1) 

For p e the matter is more delicate. Assuming that everything is smooth 
we would like to assume homogenous Dirichlet boundary condition for p £ . 
However if we consider the second equation of the system we can rewrite 
the pressure as p\ = — l/edivu 6 , from the energy priori estimate we only 
know that divu e uniformly bounded in L 2 , hence we don't have sufficient 
regularity in order to define a trace for p £ . As we will see, in the Section 
[5j we are going to decompose the pressure in two parts. One connected 
with the viscosity part of the fluid and the other one with the convective 
terms. So, our acoustic pressure waves are given by the superposition of 
two different waves with different frequencies scales. These different scales 
provide different level of regularity for p £ , therefore on the boundary we have 
to assume 

A- a p £ \ dn = 0, for any a > i. (BC2) 

This condition seems to be in contrast with the boundary data usually as- 
sociated to the pressure of the Navier Stokes equation, indeed, as it is well 
known, if one considers the Navier Stokes equation the natural behaviour 
for the pressure at the boundary is of Neumann type. Chorin in p|] dealt 
with the same issue. In p[] he considered the artificial compressibility ap- 
proximation method in the 2 — D case. Let us denote the space coordinates 
by (xi,X2) and the boundary given by the line X2 = 0, moreover denote by 

u m(i,j) = u(iAx 1 ,jAx 2 ,nAt) p^ = p(iAx 1 ,jAx 2 ,nAt) 

the approximation of the sequences u £ and p £ . In [1] Chorin assigns Dirichlet 
boundary condition to pfj in the following way, namely if the boundary line 
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X2 = is represented by j = 1, then taking into account the equation 
epf + divu £ = 0, he writes p^f 1 '- 

CP?,! 1 ~ £ Pi,i 1 = _2 Ax^^^' 2) ~~ ~ ( u l(t+l,l) ~~ U l(i-l,l))> 

where ^ are known form the previous step. 

The issue of the correct boundary condition to assign to the pressure was 
also studied by many other authors. Among them Gresho and Sani [22j . 
(see also [39], [I]) showed that Dirichlet and Neumann boundary conditions 
for the pressure give the same solution. In fact they recover appropriate 
Dirichlet boundary conditions for the pressure by taking into account the 
normal derivative of the pressure associated with the Neumann boundary 
condition and by using Green's functions. To be more precise, in our case 
we can say that in the limit the pressure plays the role of the Lagrange 
multiplier associated with the constrain div u = 0. In this regard, in analogy 
with problems of motion of constrained rigid bodies, the behaviour of the 
pressure must be deduced in terms of the velocity field u. 

Remark 3.1. For semplicity we assume homogenous boundary conditions 
(]BCip and (|BC2p . If we are in the nonhomogeneous case, let say for 
example u £ \qq = ur, by changing u £ into u £ — ur, as usual, we come into 
the case of zero Dirichlet boundary conditions. 

Now we can state our main result. The convergence of {u £ } will be 
described by analyzing the convergence of the associated Hodge decomposi- 
tion. 

Theorem 3.2. Let (u e ,p £ ) be a sequence of weak solution of the system 
(|3.ip . assume that the initial data satisfy (|IDj) and the boundary conditions 
(|BC1I) and (IBC2I) hold. Then 

(i) There exists u G L°°([0, T]; L 2 (Q)) n L 2 ([0, T]; H 1 (fi)) such that 

u £ weakly in L 2 ( [0, T] ; H 1 (O) ) . 

(ii) The gradient component Qu e of the vector field u £ satisfies 

Qu £ — ► strongly in L 2 ([0, T];L P (Q)), for any p G [4,6). 

(iii) The divergence free component Pu £ of the vector field u £ satisfies 

Pu £ — ► Pu = u strongly in L 2 ([0, T]; L 2 oc (Q)) . 

(iv) The sequence {p £ } will converge in the sense of distribution to 

p = A^div ((« • V)n) = A^Hr((Du) 2 ). 
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(v) u = Pu is a Leray weak solution to the incompressible Navier Stokes 
equation 

P(d t u - Au + (u ■ V)u) = 0, 
u(x,0) = u (x), u\ 9n = 0. 

(vi) The following energy inequality holds for all t G [0, T], 

\ [ \u(x,t)\ 2 dx + [ [ \Vu{x,t)\ 2 dxdt <i [ \u (x)\ 2 dx. (3.2) 
2 Jn Jo Jn 2 J n 

Remark 3.3. Notice that in the limit we can recover the Dirichlet boundary 
condition for u since the approximating sequence u £ G L 2 ([0, T]; -£P(0)) and 
the trace operator is continuos. 

Remark 3.4. This theorem can be easily extended to the nonhomogeneous 
equation (13. ip . by assuming 

/ e — ►/ strongly in L 2 ([0,T}; H^ 1 (ft)). 

4 Energy estimates 

In this section we wish to establish the a priori estimates, independent on 
e, for the solutions of the system (13. ip which are necessary to prove the 
Theorem 13.21 In particular we will recover the a priori estimates that come 
from the classical energy estimates related to the system (|3,1|) . 

Proposition 4.1. Let us consider the solution (u £ ,p £ ) of the Cauchy prob- 
lem for the system (|3.1|) , Assume that the hypotheses (|IDj) and the condition 
(|BCip hold, then one has 

E(t)+ [ [ \Vu £ (x,s)\ 2 dxds = E(0), (4.1) 
Jo Jn 

where we set 

m = X Qi" £ ^*)i 2 + f i^o^i 2 ) dx - ( 4 - 2 ) 



Moreover, 










is 


bounded in L°°([0, T]; L 2 (0)) ; 


(4.3) 


zp £ t 


is 


relatively compact in ff~ 1 ([0, T] x Q), 


(4.4) 


Vu £ 


is 


bounded in L 2 ([0,T] x fi), 


(4.5) 


u £ 


is 


bounded in L°°([0, T]; L 2 (^)) n L 2 ([0, T]; L 6 (0)), 


(4.6) 


(u £ -V)u £ 


is 


founded m L 2 ([0,T];L 1 (fi))nL 1 ([0,T];L 3 / 2 (fi)) I 


(4.7) 


( div u £ )u £ 


is 


founded in L 2 ([0,T\; L 1 ^)) n L\[0,T]; L 3/2 (Q)). 


(4.8) 
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Proof. We multiply, as usual, the first equation of the system (13. ip by u £ 
and the second by p £ , then we sum up and integrate by parts in space and 
time, hence we get (|4TT]) . The estimates (fO|) . (jjjjj) , (jUSj) follow from fllTTl) . 
while (|4.6p follows from (|4.ip and Sobolev embeddings theorems. Finally 
(flTTj) and gSJ) come from (|4.5f) and (|4.6|) , □ 

Remark 4.2. We want to point out as for the Navier Stokes equations, 
here, in order to get the estimate (14. ip we only used the boundary condition 
on the velocity vector field u £ , no boundary conditions are required for the 
pressure p £ . 



5 Acoustic pressure wave equation 

In order to perform our limiting process we need some more estimates that 
don't follow from the previous one. In fact from the estimates (14. 5h and 
(|4.6|) we only get the weak convergence of u £ . Here we want to exploit 
the dispersive behaviour of the pressure acoustic wave in order to damp 
its disturbing effect. Let us differentiate with respect to time the equation 
(|3.ip 2 . by using (|3.ip 1 . we get that p £ satisfies the following wave equation 

ed tt p £ - Ap £ + Adivu e - div (^(u £ ■ V) u £ + ^(divu e )u e ^ = 0. (5.1) 

Now we rescale the time variable, the velocity and the pressure in the fol- 
lowing way 

r = — , u(x, t) = u £ (x, \fer), p(x, r) = p £ (x, y/er). (5.2) 
As a consequence of this scaling the equation (|5.1|) becomes 

9 T tP — = — Adivtt + div m u ■ V) u + -(divtt)n J . (5-3) 



Now, taking into account the right hand-side of (15. 3|) we decompose p as the 
sum of two component f>\ and p2- In particular we have that p\ is related 
to the viscosity of the fluid and satisfies 

' d TT /^~ l pi - AA _1 pi = -A _1 Adivn = iq 
< A -1 pi(x,0) = <9 T A -1 pi(x,0) = 0, (5.4) 

A~ x n\aa = 0, 
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While p 2 is connected with the convective part of the fluid and verifies the 
following wave equation: 

d TT A~ 1/2 p 2 - AA~ 1/2 p 2 = A~ 1/2 div Uu • V) u + i (div u)u^j = F 2 

A- 1 / 2 p 2 (x,0) = A -1 / 2 p(x,0) d T A- 1 / 2 p 2 (x,0) = d r A _1 / 2 p(x,0), 
A-^lan = o. 

(5.5) 

Therefore we are able to prove the following theorem. 

Theorem 5.1. Let us consider the solution (u £ ,p e ) of the Cauchy problem 
for the system (|3.ip . Assume that the hypotheses (|IPp and the boundary 
conditions (IBClj) . (1BC2P /ioW. T/ien we se£ i/ie following estimate 

e3/8 \\P £ W L fw-^ + e7/8 \\ d tP £ \\ L iw^ 4 ~ V^Mhi + l|div«§||^-i 

T * t J_J 

+ || {u 6 • V)u e + i(divn e )n e |L lr 3/2 

2 iv t i>j: 



+ V / T||divu e || i 2 i a. (5.6) 



Proof. Since pi and p2 are solutions of the wave equations (|5.4|) , (|5.5|) , we can 
apply the Strichartz estimates (12. 7|) and (|2.8p . with (x,r) € x (0,T/y/F). 
First of all we use the Strichartz estimate (|2.7p with it; = A _1 pi and we get 

IIA^pilU^ + IISrA-^H^vy-M < II^iIIlilI, (5-7) 

namely 

INI^W" 2 - 4 + ll^rPi 11^4^-3,4 < div filUa^. (5.8) 

In the same way we apply the estimate (|2.8p with w = A~ x l 2 p 2 and we 
obtain 

I|A~ 1/2 P 2 ||z4 i;c + \\d T A- x l 2 p 2 \\ L%w -, A < || A- 1 / 2 ^^, 0)11^/2 + l|i ? 2|| L i 4 /2 i 

(5.9) ' 

namely, 

INI^jy-M + ||5 T p 2 || L 4w-2,4 < ||p(x,0)|U_i/2 + 119^,0)11 .-3/2 

+ \\{u- V)u+ -(divfi)fi|| rlr s/a (5.10) 
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Now by taking into account (|5.8p , (|5. 10|) it follows that p verifies 

WpWl*W- 2 ' 4 + 11^11x4^-3,4 < 11^111^4^-2,4 + \\P2\\ L 4 W -1A (5.11) 

+ \\ d rPl\\ L 4 W --iA + \\d T p 2 \\ L 4 W -2A 

< ||p(x ) 0)|U- 1/2 + ||a r p(x ) 0)|U-3/2 

n D D 



1 



+ || (u ■ V) u + -(divu)n|| x 



3/2 



lil: 

+ ^7il|divn|| L 2 L 2. (5.12) 

Using backwards the scaling (|5.2H we end up with (j5.6|) . □ 

Remark 5.2. Notice that the right hand-side of (|5.6|) is uniformly bounded 
in e by means of E(0), see Proposition 14.11 

6 Strong convergence 

This section is devoted to the proof of the strong convergence of Qu £ and 
Pu £ . In particular we will show that the gradient part of the velocity Qu £ 
converges strongly to 0, while the incompressible component of the velocity 
field Pu £ converges strongly to Pu = u, where u is the limit profile as e j 
of u £ . We start this section with some easy consequences of the a priori 
estimates established in the previous section. 

Corollary 6.1. Let us consider the solution (u £ ,p £ ) of the Cauchy problem 
for the system (|3.ip . Assume that the hypotheses (|IDp and the conditions 
(TBCip . (IBC2P hold. Then, as e | 0, one has 

ep £ — >0 strongly in L°°([0, T]; L 2 (£1)) D L 4 ([0, T]\ W~ 2 ' A (VL)), 

(6.1) 

divn £ — >0 strongly in W~ 1,oc ([0, T];L 2 (U)) D L 4 ([0, T]; W~ 3 ' 4 ({1)). 

(6.2) 

Proof. (|6.ip . (|6.2p follow from the estimates (|4.3p . (|5.6p and the second 
equation of the system (13. ip . □ 



6.1 Strong convergence of Qu £ 

Here, we wish to show that the gradient part of the velocity field Qu £ goes 
strongly to as e [ 0. This will be a consequence of the estimate (|5.6p 
provided that we observe that by using the second equation of (|3.ip we can 
rewrite Qu £ as 

Qu £ = VA^ 1 divu e = -eVA^dtp 6 . 
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Proposition 6.2. Let us consider the solution (u £ ,p £ ) of the Cauchy prob- 
lem for the system (|3.ip . Assume that the hypotheses (|ID|) and the condi- 
tions (fBCll) . (jBC2p hold. Then as e j 0, 

Qu £ — > strongly in L 2 ([0, T]; L P (Q)) for any p G [4,6) . (6.3) 



Proof. In order to prove the Proposition 16.21 we split Qu 6 as follows 

\\Q u£ \\tf t Ll < HQ« 6 -Q« e *Ja|lz^ig + l|Q« e *ia||^z5 = J l + J 2, 

where j a is the smoothing kernel defined in Lemma 12.11 Now we estimate 
separately J\ and J 2 . For J\ by using (|2.9p we get 

Ji < a 1_3 (H) ^ ||QVn £ (t)||^d^ < a~ 3 ^~v) \\QVu £ \\ L 2 L 2. (6.4) 

Hence from the identity Qu £ = —e 1 / 8 X7A~^ 1 s 7 / 8 dtp £ and by the inequality 
(12.101) we get J 2 satisfies the following estimate 

J 2 < e^VA- N ^ s d t p £ ^\\ L , LPx < e^a 2 <^)\\ey^d t p £ \\ L2tW -^ 

< e^a^-^T^y^dtp^-^. (6.5) 

Therefore, summing up (|6.4p and (16. 5p and by using (|4.5p and (|5.6p and 
remembering that Q is a bounded operator from 1? into L 2 , we conclude 
for any p G [4, 6) that 

||Qn £ || L 2 L P < Ca" 3 ^?) +C T e 1/8 a~ 2 ~ 3 ^~p). (6.6) 

Finally, we choose a in terms of e in order that the two terms in the right 
hand side of the previous inequality have the same order, namely 

a = e 1 / 18 . (6.7) 

Therefore we obtain 

\\Qu £ \\ L 2 L v < C T e^F for any p G [4,6). 

□ 

6.2 Strong convergence of Pu £ 

It remains to prove the strong compactness of the incompressible compo- 
nent of the velocity field. To achieve this goal we need to prove some time 
regularity properties of Pu £ . 
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Theorem 6.3. Let us consider the solution (u £ ,p £ ) of the Cauchy prob- 
lem for the system (|3.ip . Assume that the hypotheses (|IPP and conditions 
(jBCip . (IBC2j) hold. Then as e j 

Pu £ — > Pu, strongly in L 2 (0, T; L 2 oc (Vt)). (6.8) 

Proof. From the Proposition 14. II we know that Pu e is uniformly bounded in 
L 2 H^.. The strong convergence (16. 8p follows by applying the Theorem 12.21 
provided that for all h £ (0, 1) we have 

\\Pu e {t + h) - Pu%t)\\ L2{[0>T]xn) < C T h l /\ (6.9) 

Let us set z £ = u £ (t + h) — u e (t), then we have 



T 




\Pu £ (t + h)- Pu £ (t)\\i 2{[0iT]xn) = I I dtdx{Pz £ ) ■ (Pz £ - Pz £ *j a ) 



io Jn 

rT 



+ 



[ [ dtdx{Pz £ ) ■ (Pz e * j a ) 
Jo Jn 

h + h- (6.10) 



By using (|2.9p we can estimate I\ in the following way 



h < \\Pz £ \\ LrL 2 [ \\Pz e (t) ~ (Pz £ *j a )(t)\\ L 2dt 

Jo 

^aT^II^II^HVP^II^. (6.11) 
Let us reformulate Pz £ in integral form by using the equation (13-ip ^ hence 



h < 



[■i r rt+h 

/ dtldxl ds(Au £ - (u £ ■ V)u £ - -u e (divu e )(s,x) ■ [Pz £ * j a )(t,x) 
Jo Jn Jt 2 

(6.12) 

Then integrating by parts and by using (I2.10|) . we deduce 
h < hT l / 2 \\u £ \\ LO o L 2\\\7u £ \\ 2 2 

/ pt+h i \ V2 

+ Ca^ l2 T l l 2 \\u £ \\ LTLl [h J || (u £ • V) u £ - -(dwu £ )u £ \\ 2 Ll ds\ 

< hT l ' 2 \\u £ \\ LTLl (\Wu £ \\ Llx +Ca- 3 / 2 \\ {u £ -V)u £ - \(dwu £ )u £ \\ L 2 Ll )j 

(6.13) 

Summing up I\, I2 and by taking into account (I4,5p ~ (l4,8p . we have 

\\Pu e \t + h)- Pu £ (t)\\ 2 LH[0tT]xQ) < CT^ 2 (a + ha^ 2 + h), (6.14) 

by choosing a = h 2 ^ 5 , we end up with (16. 9p . □ 
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7 Proof of the Theorem 13.2 

(i) It follows from the estimate (|4.6D . 

(ii) It is a consequence of the Proposition 16.21 

(iii) By taking into account the decomposition u £ = Pu £ + Qu £ , the The- 
orem 16.31 and the Proposition 16.21 we have that 

Pu e _^ u strong i y in L 2 ([0) T]; Lf oc (p,)) . (7.1) 

(iv) Let us apply the Leray projector Q to the equation (|3.ip i . then it 
follows 

Vp £ = QAu £ - Q (^(u £ ■ V)u £ ) + ^u £ div Qu £ ^j . (7.2) 

Now by choosing a test function (p G Cq°(Q x [0, T]) and by taking 
into account (fl~5j) . ([63]), (pT8l) and £□]), we get, as e j 0, 

(u £ div Qu £ ,Qp) < \\Qu £ \\ L 2 L 4\\Vu £ \\ L 2 L 2\\Q^\\ LrL 4 

+ \\Qu £ \\ L 2 Li \\u £ \\ LrL 2j\VQ<p\\ L 2 Li 0, (7.3) 

(Q((u e ■ V)u £ ),tp) = {{Pu £ ■ V)Pu £ ), Qtp) + {(Qu £ ■ V)Qu £ , Qip) 

< {{Pu £ ■ V)Pu £ ,Q<p) + \\Qu £ \\ L 2 Li \\VQu £ \\ L 2 ||VQ^|| LrL 4 

- ((« • V)u), Qp) = (Q((u • V)«), y>>, (7.4) 

(QAu e , = (Pu e , AQ</>) + (Qu £ , AQy?) 

< {Pu £ ,AQip) + ||Q« £ || i?i 4||AQ^|| i?i 4/ 3 - 0, (7.5) 

where in the previous computation we took into account that Q is a 
continuos operator but doesn't commute with derivatives. So as e j 
we have, 

<Vp e , ¥>> — (VA^ 1 div((« • V)«), tp). (7.6) 

(v) In a similar way we can pass into the limit inside the system (|3.ip and 
we get u satisfies the following equation in D'([0, T] x Q) 

P{d t u - A« + (u ■ V)u) = 0, 

u(x,0) = u (x), u\ d n = 0, (7.7) 

where we used (ID) and (BC1), (BC2) and the fact that the trace 
operator is bounded (see [T7], [2]) 
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(vi) Finally we prove the energy inequality. By using the weak lower semi- 
continuity of the weak limits, the hypotheses (ID) and denoting by x 
the weak-limit of \fep £ , we have 

[ \\x\ 2 dx+ I )-\u(x,t)\ 2 dx + f [ \Vu(x,t)\ 2 dxdt 
in 2 Jq2 Jo Jn 

< liminf ( / -\u £ (x,t)\ 2 dx + [ -\p £ (x,t)\ 2 + [ [ \Vu £ (x,t)\ 2 dxdt) 
\Jn 2 Jn 2 Jo Jn J 

= liminf / \{\u £ Q {x)\ 2 + e\pl{x)\ 2 )dx= / \\u {x)\ 2 dx. 
Jn 1 Jn 1 

for any t E [0, T]. 
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